Event-by- event viscous hydrodynamics for Cu-Au collisions at s/snn = 200GeV 



Piotr Bozek 

The H. Niewodniczariski Institute of Nuclear Physics, PL-31342 Krakow, Poland 
Institute of Physics, Rzeszow University, PL-35959 Rzeszow, Poland 



Abstract 



(N 

o 

(N 



Event-by-event hydrodynamics is applied to Cu-Au collisions at ^/sjvat = 200GeV. Predictions for charged particle 
distributions in pseudorapidity, transverse momentum spectra, femtoscopy radii are given. The triangular and elliptic 
flow coefficients are calculated. The directed flow at central rapidity in the reaction plane in asymmetric collisions is 
nonzero, fluctuations of the initial profile lead to a further increase of the directed flow when measured in the event 
plane. 
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The emission of particles with soft momenta in relativis- 
i'tic heavy-ion collisions can be described in a satisfactory 
jway within the relativistic hydrodynamic model [HQ. The 
harmonic flow coefficients are found to be sensitive to the 
viscosity of the expanding fluid Fluctuations de- 

termine the asymmetry of the initial distribution in the 
i"transverse plane @, 0] and event- by- event hydrodynamics 
simulations are needed to investigate their effect on final 
spectra I3~15|- Typically collisions of symmetric nuclei 
■are studied experimentally, the most important parame- 
ters varied in these analyses are the center of mass energy, 
the size of the system and the centrality of the collision. 

An additional insight into the mechanism of particle pro- 
duction in relativistic heavy-ion collisions can be gained 
from interactions of two different nuclei, e.g. Cu-Au col- 
lisions at y/sJfN = 200GcV. The study of the multiplic- 
ity distribution and correlations in the longitudinal direc- 
tion in asymmetric collisions gives additional constraints 
jon the mechanism of the energy deposition in the early 
stage of the reaction (l6l - [2Cl |. in particular models assum- 
ing approximate boost-invariance of the initial state must 
revised. The directed flow 2lN24| in asymmetric collisions 
is partly generated by a new mechanism, as discussed lat- 
ter in this Letter. The measurement of differences between 
the directed flow of proton and antiprotons could give in- 
formations on the different rate of baryon stopping on the 
Cu and Au side of the fireball which could be compared 
to model predictions 251 27t. The reaction plane corre- 
lations in the fireball |28l l29j are modified in asymmet- 
ric collisions, with stronger correlations between odd and 
even order event planes. As the size and the energy den- 
sity in the fireball are sufficient to allow for a substantial 
phase of collective expansion, we apply the event- by- event 
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Figure 1: The eccentricity 62 (solid line), triangularity £3 (dotted 
line) and the dipole asymmetry t\ (dashed line) of the initial fireball 
as function of the number of participant nucleons in Cu-Au collisions 
from the Glauber Monte Carlo model |3d(l . The vertical lines separate 
different centrality classes. 



viscous hydrodynamic model to obtain predictions for ba- 
sic observables for particles of soft momenta produced in 
asymmetric Cu-Au interactions. 

The second order viscous hydrodynamic equations are 
solved in 3 + 1-dimensions [IH [3l[ . We use two values for 
the shear viscosity coefficient r\j s = 0.08 or 0.f6 (unless 
stated otherwise the default value is 0.08), the bulk vis- 
cosity is (/s — 0.04 in the hadronic phase. At freeze-out 
particle emission is performed using a Monte Carlo proce- 
dure |32] , with a freeze-out temperature of 150MeV. This 
freeze-out temperature gives a satisfactory description of 
particle spectra in Au-Au collisions 12J. The hydrody- 
namic expansion is followed in each event with random 
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initial density. The initial density is generated from the 
Glauber Monte Carlo model [3(| with a Gaussian wound- 
ing profile for nucleon-nucleon collisions [33} . The entropy 
density in the transverse plane x-y and space-time rapidity 
r)\\ is given as a sum over the participant nucleons 



weighted with Gaussians gi(x,y) of width 0.4fm centered 
at the positions of the participant nucleons i, a term from 
binary collisions with a = 0.125 is included, Nf oU is the 
number of collisions of the participant nucleon i. The form 
of the longitudinal profile f± and the factor d = 2.5 GeV 
are taken the same as for Au-Au collisions at 200GeV 



In Fig. [T]are shown the initial eccentricity £2, triangularity 
€3 and dipole asymmetry t\ of the fireball as function of the 
number of participants, averaged over the initial density 
with weights r 2 , r 3 and r 3 respectively. In the following we 
present results for two centrality classes 0-5% and 20-30% 
defined by the number of participants 190 < N part and 
81 < Np ar t < 114 respectively. For each centrality class 
considered we generate 100 hydrodynamic events, and for 
each hydrodynamically generated freeze-out configuration 
1000-1500 statistical emission events. The relative statisti- 
cal error on V2 and V3 from sampling the initial eccentricity 
and triangularity distributions with 100 events is 4 — 5%. 
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Figure 2: Distribution of charged particles in pseudorapidity for 0- 
5% (solid line) and 20-30% (dashed line) centrality classes. The Au 
momentum is directed towards positive rapidity. 

The distribution of charged particles in pseudorapidity 
dNch/dr/ps is asymmetric (Fig. [2]). It reflects the asym- 
metry of the initial entropy density (Eq. [TJ [ll, Q • The 
number of participant nucleons is larger in the Au nucleus, 
with the asymmetry increasing for central collisions. Using 
the initial condition extrapolated from Au-Au interactions 
at the same energy, the predicted multiplicities at central 
rapidity are dN c h/dr]ps — 320 and 140 in 0-5% and 20- 
30% centrality classes respectively. 
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Figure 3: Transverse momentum spectra of 7r+ (solid lines), K + 
(dashed lines) and protons (dotted lines) for 0-5% and 20-30% cen- 
trality classes (upper and lower curves respectively). 



The transverse momentum spectra of ir + , K + and pro- 
tons shown in Fig. [3] are softer than in Au-Au collisions. 
The predicted average transverse momenta of pions, kaons 
and protons in central collisions are 405, 587 and 775McV 
respectively. To obtain correct particle ratios nonzero 
chemical potentials are introduced at freeze-out, although 
the equation of state used to calculated the dynamics of 
the system is at zero baryon density. The equation of state 
is moderately changing with /1 at small baryon densities 
[3"r3 | , and the equation of state at zero baryon density can 
be used for central rapidities in collisions at 200GeV. The 
the energy in the Cooper-Frye formula at freeze-out is con- 
served to the order /i 2 /T 2 . The baryon chemical potential 
fiB = 22MeV at the freeze-out temperature 150MeV as- 
sures the same ratio p/p as measured for Au-Au interac- 
tions [35} . 

The elliptic and triangular flow coefficients as function 
of p± are presented in Fig. 0] The differential flow co- 
efficients are calculated using the two-particle cumulant 
method [37[. The elliptic flow coefficient for semi-central 
events is significantly larger than in central events. It re- 
flects the larger value of the initial eccentricity in the 20- 
30% centrality class (Fig. [IJ. The triangular flow is sim- 
ilar in the two centrality classes studied. With increasing 
value of shear viscosity the ellip tic and triangular asymme- 
try coefficients are reduced 13J, |38j . The influence of shear 



viscosity is larger in peripheral events and for the triangu- 
lar flow i>3 3^, 4(j • The integrated value of the elliptic flow 
of charged particles in the p± range 150-2000MeV is 0.022 
(0.020) for central and 0.048 (0.043) for peripheral events 
when 77/s = 0.08 (0.16). The corresponding v 3 coefficients 
for charged particles are 0.013 (0.011) and 0.016 (0.012). 
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Figure 4: Elliptic (panel a) and triangular (panel b) flow coefficients 
of charged particles as function of transverse momentum for cen- 
tralities 0-5% (solid lines) and 20-30% (dashed lines). The results 
of the calculation with rj/s = 0.08 are presented with lines and for 
•q/s = 0.16 using lines with triangles. 



In collisions of symmetric nuclei at ultrarelativistic en- 
ergies the directed flow exhibits a component odd in pseu- 
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dorapidity [21|, |22|, |41j, |42j ■ At central rapidities nonzero 
directed flow occurs due to event-by-event fluctuations of 
the density profile 11, 43-45}. In the fireball created in 
the collision of two asymmetric nuclei such as Cu-Au col- 
lisions the profile is asymmetric in the reaction plane also 
on average. When defining the reaction plane always with 
the Au nucleus in the x > half plane a nonzero compo- 
nent of the directed flow with respect to the reaction plane 
appears. 

In Fig. [5] is shown the px integrated directed flow coef- 
ficient v\ of charged particles with respect to the reaction 
plane as function of pseudorapidity for centrality 20-30%. 
At central rapidities the directed flow is negative, more 
particles flow in the direction of the Cu half plane. Be- 
sides the nonzero component even in rjps a smaller odd 
component of v\ is visible in Fig. [5] In the hydrodynamic 
model the odd component of the directed flow is generated 
from the expansion of the fireball tilted away from the col- 
lision axis [2l], H^. The odd component of the directed 
flow is significantly reduced when viscosity increases, this 
effect comes from the reduction of the longitudinal pres- 
sure in the early phase of the collisions for higher viscosity 
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Figure 5: Directed flow of charged particles as function of pseu- 
dorapidity with respect to the reaction plane. The results of the 
calculation with -q/s = 0.08 are presented using a solid line and for 
rj/s = 0.16 using a solid line with triangles. 



The even component of the directed flow as function of 
transverse momentum in the reaction plane 



Vl ( P± ){RP} = (cos(& - *bp)) 



(2) 



is estimated as the average for charged particles with 
\vps\ < 1- The direction of the reaction plane ^pp is 
chosen in the direction of the Au nucleus. Defining the 
reaction plane using all the nucleons in the Au and Cu 
nuclei or only the spectator nucleons in the two nuclei in 
each event gives indistinguishable results. The directed 
flow v\(p±) is negative for small momenta and changes 
sign for p± ~ 850MeV (Fig. [B]). The component of the 
directed flow odd in pseudorapidity 

wi(p x ){i?P}(odd) = (sgn(r]ps) cos(& - V RP )) (3) 



is much smaller than the even one (we take charged parti- 
cles with \t]ps\ < 2). 

Fluctuations of the fireball density in each event change 
the orientation and the magnitude of the directed flow in 
each event 



431 ] . We follow the procedure of Refs. [11|.|47| 



where the Q vector of the event plane is defined with a 
weight reducing the contribution of momentum conserva- 
tion to the directed flow 



with Wi = p± - (p 2 L )/(p±). 
directed flow coefficient is 



(4) 

The Q weighted value of the 



v^p^EP} = 



(Qcos((fc - fri)) 



(5) 



The directed flow coefficient vi(p±) with respect to the 
event plane has the same form as the even component 
defined in the reaction plane, but with a slightly larger 
magnitude (solid line in Fig. [BJ). This means that fluctu- 
ations increase the directed flow at central rapidity. The 
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Figure 6: Directed flow of charged particles as function of transverse 
momentum with respect to the reaction plane (dashed line) and the 
event plane (solid lines) for |»?ps < 1, the dots represent the directed 
flow with respect to the reaction plane odd in pseudorapidity for 
|?7ps| < 2. The results of the calculation of v\{EP} with r]/s = 0.16 
are presented using a solid line with triangles. 
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calculated directed flow in the event plane does not depend 
strongly on the value of shear viscosity. 

The correlation function for same charge pion pairs is 
calculated using the momenta and positions of pions emit- 
ted in each event [HI, |48[ . From the correlation function 
in relative momentum of the pair the femtoscopy radii are 
extracted. In Fig. 0is shown the dependence of the three 
radii R ou t, R S ide, and Ri on g on the average pion momen- 
tum. This dependence is similar as seen in Au-Au colli- 
sions, but with smaller values of the radii. In addition, 
we calculate the femtoscopy radii using one hydrodynamic 
simulation starting from the average initial condition cor- 
responding to centrality 0-5% (solid lines in Fig [7]). The 
results are very similar as obtained in event by event sim- 
ulations. It shows that flow fluctuations in the event by 
event evolution are too small to affect significantly the 
femtoscopy radii. 

We present predictions of the event by event viscous 
hydrodynamic model for Cu-Au collisions at y / 5jv r /v = 
200GeV. The charged particle multiplicity and femtoscopy 
radii are smaller than in Au-Au interactions, reflecting the 
reduced size of the system. The transverse momentum 
spectra are steeper and the elliptic and triangular flow 
smaller than for the Au-Au system. 

A novel aspect for collisions of asymmetric nuclei is seen 
in the directed flow. The asymmetry of the fireball density 
in the reaction plane leads to the formation of a directed 
flow of charged particles at central rapidities. This contri- 
bution appears when the orientation of the pair of Cu and 
Au nuclei in the reaction plane is determined in each event 
and is nonzero also in the absence of density fluctuations. 
Density fluctuations further increase the directed flow at 
central rapidities when measured in the reaction plane. 
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Figure 7: The femtoscopy radii R ou t, R s ide> an d Riong as function 
of the average transverse momentum of the pion pair. The dashed 
lines represent the results of the event-by-event hydrodynamics and 
the solid lines the results obtained in a simulation using one average 
initial condition. 
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